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ABSTRACT 

Radio relics in galaxy clusters are associated with powerful shocks that (re)accelerate rel¬ 
ativistic electrons. It is widely believed that the acceleration proceeds via diffusive shock 
acceleration. In the framework of thermal leakage, the ratio of the energy in relativistic elec¬ 
trons to the energy in relativistic protons should should be smaller than K e / p ~ 10 -2 . The 
relativistic protons interact with the thermal gas to produce 7-rays in hadronic interactions. 
Combining observations of radio relics with upper limits from 7-ray observatories can con¬ 
strain the ratio K e / p . In this work we selected 10 galaxy clusters that contain double radio 
relics, and derive new upper limits from the stacking of 7-ray observations by Fermi. We 
modelled the propagation of shocks using a semi-analytical model, where we assumed a sim¬ 
ple geometry for shocks and that cosmic ray protons are trapped in the intracluster medium. 
Our analysis shows that diffusive shock acceleration has difficulties in matching simultane¬ 
ously the observed radio emission and the constraints imposed by Fermi, unless the magnetic 
field in relics is unrealistically large (2> 10 //G). In all investigated cases (also including re¬ 
alistic variations of our basic model and the effect of re-acceleration) the mean emission of 
the sample is of the order of the stacking limit by Fermi, or larger. These findings put tension 
on the commonly adopted model for the powering of radio relics, and imply that the relative 
acceleration efficiency of electrons and protons is at odds with predictions of diffusive shock 
acceleration, requiring K e / p ^ 10 — 10~ 2 . 

Key words: Galaxy clusters; intergalactic medium; shock waves; acceleration of particles; 
gamma-rays. 


1 INTRODUCTION 

Radio relics are steep-spectrum radio sources that are usually de¬ 
tected in the outer parts of galaxy clusters, at distances of ~ 0.5 — 3 
Mpc from their centres. They are very often found in clusters 
with a perturbed dynamical state. There is good evidence for their 
asso ciation wit h po werful merger shocks, as early suggested by 
lEnsslin et alj d 1998h . Among them are giant double-relics that 
show t wo large sources on opposite sides of the host cluster’s centre 
(e.g., iFeretti et ahll 2012 l : iBonafede et alJ I20T2I : Ide Gasperin et all 
120141) . These relics are associated with shocks in the intracluster 
medium that occur in the course of cluster mergers. Only a tiny 
fraction of the kinetic power dissipated by typical cluster merger 
shocks (<C 10 3 ) is necess ary to po we r the relics, and di ffu¬ 
sive s hock acceleration (DSA, Canriolil e.g. |20 1 2l : Kang & Ryu] e.g. 
120131 for modern reviews) has so far been singled out as the most 
likely mechanism to produce the relativistic electrons and to pro- 
duce the observed power-law radio spectra (e.g. iHoeft & Brtiggenl 
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200 7). However, if the standard DSA model is correct, the same 
process should also lead to the acceleration of cosmic-ray (CR) pro¬ 
tons. Indeed, the process should be much more efficient for protons, 
owing to their larger Larmor radius Q. 

To date, high-energy observations of nearby galaxy clusters have 
not revealed any diffuse 7 -ray emission resulting from the interac¬ 
tion between relativisti c protons and thermal particles of the intra- 
cluster mediu m (ICM ) 1 Reimer et alJl2003l: Aharoni^etalJl^OOSl; 
Aleksic et alj 1201(1 12012 : Arlen et al. 2012 : Zandanel & Andol 


20141) . Recently, the non-detection of diffuse 7 -ray emission from 


clusters by Fermi has put the lowest upper limits on the density of 
CRs in the ICM, < a fe w percent of the thermal gas en ergy within 
the clusters virial radius dAckermann et al.l t2010. 2014). Moreover, 
the stacking of subsets of cluster observations leads to even lower 


1 Instead, since the Larmor radius of thermal electrons is much smaller 
than the typical shock thickness, thermal electrons cannot be easily accel¬ 
erated to relativistic energies by DSA. T his so c alled i njecti on problem for 
electrons is still largely unresolved fe.g. lBrunetti & Jonesi20l4) 
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upper-limits (Huber et alj2013l :lGriffin et al J2014h . These low lim¬ 
its on the energy content of CR s can be used to constrain shock- 
acceleration models. Recently, Ivazza & Briiggenl C20I 4) have al¬ 
ready suggested that the present statistics of radio observations, 
combined with available upper limits by Fermi places constraints 
on PS A as the source of giant radio relics. In Vazza _&Briiggen 
< hoi4 we assumed that the population of clusters with radio relics 
was similar to the population of (non-cool-core) clusters for which 
the stacking of Fermi clusters was available. In the present paper, 
we repeat a similar analysis by comparing to a more realistic stack¬ 
ing of the Fermi data. Our method is outlined in Sec. 12.II while our 
results are given in Sec. [3] In the latter Section, we also discuss on 
the role played by the several open parameters in our modeling. We 
find (Sec. [3} that the present upper limits from Fermi imply energy 
densities of CR-protons that are too low to be explained by stan¬ 
dard DSA: if DSA produces the electrons in relics, then we should 
have already detected hadronic 7 -ray emission in some clusters, or 
in stacked samples. In our conclusions (Sec. 0. we discuss pos¬ 
sible solutions to this problem, as suggested by recent hybrid and 
particle-in-cell simulations of weak, collisionless shocks. 


2 METHODS 


2.1 Semi-analytical cluster mergers 

Our aim is to test DSA by making quantitative predictions of 
the hadronic 7 -ray flux assuming that the CR-protons come 
from the same shocks that produce radio relics (Sec. 12.2. H . 
Semi-analytical methods with initial conditions tuned to match 
observable param eters of radio halos have been widely used in 
the literature (e.g. lGabici & Blasi 20031: ICassano & Brunettill2005l : 


iHoeft & Briiggenl 20071) . They have obvious limitations owing 
to the lack of 3-dimensional detail and the crude geometrical 
assumption on the merger scenario (i.e. spherical symmetry and 
simple radial distribution for the gas). Still, for this specific 
problem they are helpful, as the energy density of CR-protons 
should be a simple function of the shock parameters and of the 
volume crossed by each shock wave. 


Our approach is similar to the method used in 
IVazza & Briiggenl J2014I) : for each cluster in our sample, we 
model the shock trajectories and the associated CR-acceleration. 
We use simple ID models of shock propagation in stratified 
atmospheres and use observables from radio and X-ray data, such 
as the spectral index, the distance from the centre, the radio power 
and the largest linear scale of e ach relic iBonafede et al.1 d2012l) 
(see also lde Gasperin e t alJ 201 41), Below we s ummarise the most 
important steps from Vazza & Briigge^ d2014l) (see also Fig. 0, 
while in Sec l3.2l we discuss the most relevant uncertainties in our 
model and assess their impact on our results: 

(i) We infer the Mach number, At, from the spectral index of 
the radio spectrum, s, at the injection region via At = </ |±| (J is 
the slope of the particle energy spectra and S = 2s). This assumes 
that the radio spectrum is dominated by the freshly injected CR- 
electrons at the shock, an assumption that might be poor for spectra 
integrated over large downstream volumes (in which cases the radio 
spectrum is steeper by 0.5). Bootstrapping with random deviates 
from the Mach number thus determined can quantify the errors (see 
Sec.0. 

(ii) The upstream (i.e. pre-shock) gas density, n u , at the relic 


is computed using a , 0 -model profile for each host cluster, with 
P = 0.75 and the core radii scaling as r c = r c ,Coma(Tcoma/T , ) 1/2 
(which follows from the self-similar scaling), where r c ,Coma = 
290 kpc. This is the only way to regularise our dataset, as a more 
detailed reconstruction of the gas density from the literature is 
missing for most of these objects. Both are clearly a simplifica¬ 
tion but in reality we probably find higher gas densities and tem¬ 
peratures along the merger axis, due to clumping and enhanced gas 
compression which will only increase the 7 -ray emission. More¬ 
over, the propagation of each merger shock is determined by the 
temperature in front of it, i.e. in the upstream region, which present 
X-ray observation can hardly constrain for any of these objects. 
Therefore we always consider an upstream gas temperature, T u , 
based on the Lx — T500 relation for each host cluster dPratt et alJ 
120091) . 


Kang & Ryu 


Kang & Rvu 


(iii) We compute the kinetic power for each shock as <Fkin = 
n u VsS/2, where v s = Alc s (c s oc \EE). 

(iv) We assume that a fraction of the kinetic power goes into 

CR-protons <f>cRp = T](M)i E>kin, where the efficiency, is 

a non-linear function of the Mach number. It has been derived 
for several DSA scenar ios (e.g. Kang & Jones (20071 : 
l2013| jHong et alJl2014l) . Here we use r/(M) given by ] 

00131) . which were estimated based on simulations of nonlinear 
DSA, considering an upstream P = 100 plasma and including a 
phenomenological model for the magnetic field amplification and 
Alfvenic drift in the shock precursor, due to accelerated CRs. This 
function predicts an acceleration efficiency of ~ 1 percent for 
At =_3, steeply rising to ~ 10 percent for At = 5 (see Fig. 4 
of iKang & Rvul20l3l) . The corresponding power into CR-electrons 
is set by assuming a fixed electron-to-proton ratio, K e / p = 0.01: 
dcRe = A' e / p t)(M)<I>kin. This ratio is already conservative, as re¬ 
cent models of particle acceleration in supernovae suggest an even 
lower value, K e / p ~ 10 -3 dPark et akll2014l) . which would result 
into a 10 times larger hadronic 7 -ray flux than for our value of 


Ke/p- 

(v) The magnetic field at the r elic, B, is derived from t he radio 
power via the equations given by iHoeft & Briiggeiil d2007t) . In this 
model the monochromatic radio power at frequency u, P„, depends 
on the shock surface area, S (which we derive from the projected 
size of the relic, assuming that the relic has a circular shape), the 
upstream electron density, n e (computed from n u by assuming a 
mean molecular mass of 0 . 6 ), the upstream electron temperature, 
T e « T u , the spectral index of the radio emission, s, and the relic 
magnetic field, B re Uc, in the following way: 


dP 6.4 • 10 34 erg 


dv 


s ■ Hz 


• S-n e -v(M)K e/F 


B 




,.s/2 f>2 1 d 2 

V D CMB + -“relic 


,( 1 ) 


where Bomb is the equivalent field of the cosmic microwave back- 
ground@In our simplest model ("basic model") we let the magnetic 
field as a free parameter without upper limits, while in a more real¬ 
istic scenario ("Bcap” scenario) we imposed a maximum magnetic 
field of Breiic.max = 10 [iG for all relics (Sec l3.2.j~l l. 

(vi) Explaining the observed radio emission from At ^ 3 
shocks is a problem, as the required electron acceleration ef¬ 
ficiency at these weak shocks can become unrealistically large 
GMacario et all DPI II) . It has been suggested that the contribution 
from shock re-accelerated electrons can alleviate this problem, as 


2 We notice that, compared to the original Equation by IHoeft & Briiggenl 
<2007t) . we use here upstream values for density and temperature, as the 
function rj(M) also accounts for the shock compression factors. 
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it would mimic the effect of having a h igher acceleration efficiency 
d Kan si et al. 2012j: Pinzke et alii20131) . To model re-acceleration, 
as in IVazza & Brtigeenl ll20 1 4) . we used an increased CR-proton 
acceleration effici ency as a function of Ma ch number, following 
iKang et aT] d2007l) and iKang & Rvul d2013l ) who showed that the 
net effect of re-accelerated and freshly injected cosmic rays can be 
modeled by a rescaled efficiency r/(M). This depends on the energy 
ratio between pre-existing cosmic rays (-Ecr) and the thermal gas 
(E g ). In the following we parametrize this ratio using the param¬ 
eter e = Ecu/Eg, and explore the cases e = 0 (single injection, 
no-reaccelerated electrons), and bracket the trend of re-acceleration 
using e = 0.01 and t = 0.05. Notice that the latest limits from 
Fermi only allow s e ^ a few percent, for flat radial distributions of 
CRs dAckermann et alJl2014l) . 

(vii) In order to compute the spectrum of electr ons i n the case 
of a re-accelerated pool of CRs, we follow [Kang et all d2012l ) 
and assume that the blend of several populations of pre-existing 
CR-electrons are characterised by a power law with index 5 e = 
(4<5 + l)/3, where S is the spectral index of the energy spectrum de¬ 
rived at the relic as before. As in lKang et al.l d2012t) we fix the spec¬ 
trum of re-accelerated CRs to the particle spectrum at the relic in 
those cases where the spectrum is flatter than the (steep) spectrum 
associated with the M = 2 shock that is supposed to re-accelerate 
them. 

(viii) Once the shock parameters are fixed, we estimate the en¬ 
ergy injected in CR-protons that are assumed to stay where they 
have been predicted. We assume that each shock surface scales with 
the cluster-centric distance, r, as S(r) = So(r / R le n c ) 2 (i.e. the 
lateral extent of the shock surface is set to largest linear size of the 
relic, which decreases with oc r inwards). The cumulative energy 
dissipated into CR-protons is given by 


C^relic 


E CR = 


'I’crpMvs dr = 


^relic 


ri(M)n u (r)Vs S(r) 


dr.(2) 


The shock surface and strength will vary with radius, and so will 
the CR-proton acceleration efficiency, r/(M). Our fiducial model 
assumes that the Mach number of shocks released by the merger is 
constant across the whole volume of interest, while in Sec. l3.2.3l we 
also test a scenario in which the Mach number scales with radius. 
The lower integration limit in the equation for Ecr is the core ra¬ 
dius, r c , meaning that the shocks are assumed to be launched only 
outside of the clus t er core, which is sup ported by simulations (e.g. 
IVazza et alj| 2012 fcl : lskillman et al .1120 1 3l) 

(ix) We compute the hadronic 7 -ray emiss ion, J 7 , following 
dPfrommer & Enfili5l2004l : iDonnert et al.ll2010lf .with the only dif¬ 
ference that for the hadronic cross-section we use t he parametri- 
sation of the proton-proton cro ss section given by iKelner et alJ 
l2006h . as in iHuber et~al] J20 1 3l ). In detail, we compute for each 
radius the source function of 7 -rays as: 


Q-yiE-y) 


2 4 ^ a pp (E)cn e (r)K P /r , Ain 

3<5 7 <5-1 1 j GeV 


m^oc 2 
GeV ’ 


(3) 


where n e is the upstream electron density, computed from the 
gas density by assuming a molecular mean weight /x = 0.6. The 
spectrum of the 7 -ray emission depends on the assumed Mach 
number across the cluster, and therefore is either a function of 
M(r) in the single acceleration model, or a constant in the re¬ 
acceleration case. Once the spectral index, S, of the particle spec¬ 
trum is fixed, the spe ctrum of the 7 -ray emiss ion is given by 
<5 7 = 4(5 - 1/2)/3 liPfrommer fc Enfilinl 12004 ). Here, we con¬ 
sider hadronic emission in the energy range [0.2-300] GeV, which 
are compared to the stacked emission from Fermi described in the 


following section (Sec. 12.2.21 . m n and m p are the masses of the n° 
and the proton, respectively. The threshold proton energy is taken 
as Emi n = 780 MeV and the maximum is i? max = 10° MeV (the 
actual value is actually irrelevant, given the steep spectra of our 
objects) The effective cross-section we used is 


a pp (E) 


(34.3+ 1.88L + 0.25L 2 ) 



(4) 


withL = \n(E/\ TeV) wdEth = m p + 2m w + mj/2m p ~ 1.22 
GeV (Eq. 79 of iKelner et all20061) . The normalisation factor K p 
is: 


K p = 


(2 — 5)Ecn{r) 

( rp—8+2 771—<5+2\ * 

V^max ■ C/ min ) 


(5) 


The emission per unit of volume, A 7 , is obtained by integrating 
the source function over the energy range: 


E 2 

A 7 (r) = f dE 7 q 7 (E 7 ) 

El 

_ a-ppm^c 3 n u (r)K p (E min )~ S E„,;„ 




/ i 7 + l 5 7 — 1 
l 2 i5 7 ’ 2 S 7 




( 6 ) 


where B x (a, b) denotes the incomplete /3-function and [/(a;)]/ = 
f(a) — f(b). Integrations over radii are performed by summing 
over radial shells of thickness 10 kpc. Finally, the hadronic 7 -ray 
emission for each radial shell along the trajectory of the shock 
is given by A 7 (r)S(r) dr and the total hadronic emission in the 
downstream region of each relic is given by the integral, I 7 = 
I*™' 1 ' A fr)S(r) dr. 

The above set of approximations minimises the hadronic 7 -ray 
emission from clusters: the presence of gas clumping and substruc¬ 
ture is neglected, the injection of CRs is limited to regions outside 
of the dense cluster cores, and additional acceleration of CRs by 
earlier shocks, turbulence, supernovae and AGN is not taken into 
account. A number of assumptions have to be made in the previous 
steps. We discuss all most important in Sec. 13.21 and explore the 
effects uncertainties in the parameters of the model. 


2.2 Observations 

2.2.1 Radio data 


We restrict our analysis of to double radio relics, as these sys- 
tems are caused mo s t clearly by major merger events (e.g. 


terns are caused mo s t clearly by m aior merger event s (e.g. 
Roettieer et alJ 1 19991 : Ivan Weeren et all 1201 lal : Iskillman et alJ 


20 1 3l) . and they should be less affected by projection effects be¬ 


cause of large cluster-centric distances dVazza et al.l120 1 2 bl). We 
select double-relic sources from the collection of Bonafede et al. 
J2012I) and further restrict analysis to the sources at high Galactic 
latitude (|fe| > 15°) to avoid strong contamination by the bright dif¬ 
fuse 7 -ray emission from the Galactic plane (see below). The final 
sample is made of 20 relics from 10 clusters: MACSJ1752, A3667, 
A3376, A1240, A2345, A3365, MACS J1149, MACSJ1752, 


PLCKG287, ZwClJ2341 and RXCJ1314. The values of radio pa¬ 
rameters for these objects (e.g. total power, radio spectral slope, s, 
largest linear scale of each object and distance from the centre of 
the host cluster) are given in Tab. 1. 
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Figure 1. Schematic view of our method for computing 7 -ray emission 
from accelerated CR-protons downstream of the double relics. 


2.2.2 7 -ray data 

We analysed 7 -ray data collected by the Large Area Telescope 
(LAT) on board Fermi (hereafter Fermi-LAZT). Fermi-LAT 
lAtwood et al .2009h is a pair-conversion 7 -ray telescope operating 
in the 20 MeV - 300 GeV band. We collected all the data obtained 
during the period 2008-08-04 to 2013-11-01 and analysed them 
using the Fermi Science Tools software package v9r32p5 and the 
P7SOURCE_V6 instrument response files. For each source listed 
in Table 1, we extracted a rectangular region of interest (ROI) with 
side ~ 20° x 20° centred on the source. For each source, we then 
constructed a model including the Galactic emission^ the isotropic 
diffuse background Q and all sources listed in the 2FGL catalog 
in a radius of 20 ° around our source, with spectral models and 
parameters fixed to the values given in 2FGL. We then performed a 
binned likelihood analysis dMattox et aljl 199bll on each individual 
ROI to determine the parameters of the 7 -ray emission model. The 
normalization of each component (diffuse background components 
and point sources) was left free to vary during the fitting procedure. 

As an alternative hypothesis, we then added to the model a 
pointlike test source fixed to the cluster position. The test source 
had a spectrum computed from the simulated CR proton popu¬ 
lation, accelerated at shocks (Sec. 12.H . which should be injected 
by M ~ 3 shocks in most cases. This spectrum i£ the best 
g uess from our cosmological numerical simulations dVazza et~al] 
2014b ) and is close t o the ’’universal CR-spectrum” suggested by 
Pinzke & Pfrommerl d2010l) . The resulting photon spectrum was 
discretised and used as a template for the expected emission of the 
shock-accelerated CR protons. This spectral model closely resem¬ 
bles a power law with a photon index of ~ 2.6 at energies > 1 GeV 
and is significantly flatter at lower energies because of the strong 


3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/aux/gal_2yearp7v6_v0.fits 

4 http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/aux/iso_p7v6source.txt 



Figure 2. Co-added Fermi-LAT count map in the 200 MeV - 300 GeV 
energy range for all the clusters listed in Table 12.2.21 The circles indicate 
apertures of 2.5° and 5° around the stacked cluster position to highlight the 
source and background regions, respectively. 


dependence of the p — p interaction cross section near the thresh¬ 
old energy for pion production (see Appendix A. 1 of iHuber et~aD 
120131) . We then estimated the test statistic (TS) defined as 


TS = -2 In 


£ max 
0 

Qrnax ’ 


(7) 


where C™ ax and are the maximum likelihood values for 

the null and alternative hypotheses, respectively. No significant 
signal (TS > 25) was observed for any of the clusters from Table 
1. Therefore, we computed upper limi ts at the 95% confidence 
level (CL) for all sources as in lAckermann et all J2014I) . which we 
report in Table 3 for two energy ranges (0.2 — 300 and 1 — 300 
GeV). The choice of a pointlike source at the cluster positions is 
not critical. The mean spectra assumed for the relics are steep, and 
most of the photons should come from the low energy end of the 
spectrum, where the resolution of Fermi is too coarse to distinguish 
between poinlike and extended sources in our objects (~ 4° at 200 
MeV). 


In order to lower these limits, we performed a stacking anal- 
ys is of the sou r ces in our sa mple following the method presented 
in IHuber et al.l d 2012 L I20I3I) . The stacking of the sources is per¬ 
formed by adding step by step the individual ROIs after having 
simulated and subtracted the surrounding point sources. The result¬ 
ing co-added map for our stacked sample is given in Fig. [2] The co¬ 
added data are then fit using the same source+background model as 
described above. Again, no significant emission is measured within 
the stacked volume of these clusters (T5 a tacked = 0 . 21 ), which 
results in an upper limit of T’mean < 1-90 • 10 ~ 10 ph/cm 2 /s (95% 
CL, 0.2 — 300 GeV band) to the average flux per cluster. For more 
details on the stacking procedure and a thorough validation of the 
method using simulated data, we refer the reader to IHuber et alj 

d2012h . 
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Table 1. Main observational parameters for the radio relics and clusters considered in this paper: redshift (2nd column). X-ray luminosity (3rd), distance from 
the cluster centre for each relic in the pair 4th and 5th column), largest linear scale (6-7), radio power (8-9), radio spectral index (10-11) of relics, u pper limits 
of 7 -ra y e mission at 0.2 — 300 an d 1 — 300 GeV for each host cluster (12-13) and value of the test statistics (Eq.|7). The radio data are taken from Fer etti et al. 
(m! and lBonafede et all J2012I) , while the 7 -ray data have been derived from the Fermi catalogue in this work. 


object 

z 

Lx 

10 44 erg/s 

ri 

Mpc 

Mpc 

Ri 

Mpc 

7?2 

Mpc 

lo gio( p R,i) 

[erg/s] 

lo Slo( P R,l) 

[erg/s] 

Si 

S2 

logio(GL 0 .2-30o) 

[ph/(scm 2 )] 

i°gio( pp i-30o) 

[ph/(scm 2 )] 

TS 

A3376 

0.047 

1.04 

0.80 

0.95 

1.43 

0.52 

40.026 

39.936 

1.20 

1.20 

-8.65 

-9.71 

1.36 

A3365 

0.093 

0.41 

0.56 

0.23 

1.00 

0.70 

40.086 

39.146 

1.55 

1.93 

-8.61 

-9.66 

6.17 

A1240 

0.159 

0.48 

0.64 

1.25 

0.70 

1.10 

39.748 

39.991 

1.20 

1.30 

-9.21 

-9.89 

0 

A2345 

0.176 

2.56 

1.50 

1.15 

0.89 

1.00 

40.544 

40.577 

1.30 

1.50 

-8.55 

-9.55 

6.19 

RXCJ1314 

0.244 

5.26 

0.91 

0.91 

0.57 

0.94 

40.350 

40.714 

1.40 

1.40 

-9.06 

-9.94 

0 

MACSJ1149 

0.540 

6.75 

0.82 

0.76 

1.39 

1.14 

40.879 

41.100 

1.20 

1.42 

-9.53 

-10.11 

0 

MACSJ1752 

0.366 

3.95 

1.34 

0.86 

1.13 

0.91 

41.651 

41.292 

1.21 

1.13 

-8.92 

-9.98 

0.04 

A3667 

0.046 

1.02 

1.42 

1.83 

1.30 

1.30 

40.044 

39.945 

1.39 

1.31 

-8.59 

-9.64 

3.16 

ZwCLJ2341 

0.270 

2.00 

0.25 

1.20 

1.18 

0.76 

40.726 

40.377 

1.90 

1.22 

-9.41 

-10.35 

0 

PLCKG287 

0.390 

8.29 

1.93 

1.62 

1.58 

3.00 

41.401 

41.322 

1.26 

1.54 

-8.81 

-9.68 

0.49 





Figure 3. ID profiles of various quantities along the projected propagation 
radius of the two relics in A3667, as assumed or predicted by the basic 
model (Sec. 12 TV First panel: gas density profile. Second panel: profile of 
the cumulative shock energy and CRp-energy dissipated in the downstream. 
The additional numbers in colours give the magnetic field necessary to re¬ 
produce the observed radio power using Eq.[l] Third panel: predicted 7 -ray 
emission. The vertical lines delimit the regions where the shocks have been 
launched, while the horizontal line gives the single-object limit from Fermi. 


3 RESULTS 
3.1 Basic model 

A typical set of predictions from our baseline model is given in 
Fig. a for the two relics in A3667. The integrated kinetic energy 
that crosses the shock (upper lines) and dissipated CR-proton 
energy (lower lines) in the second panel are integrated along the 


shock trajectory. In the single injection model (e = 0) the predicted 
acceleration of CRs is very different for the two relics, and follows 
from the different assumed Mach numbers: M = 3.8 for the 
northern relic and M = 1.7 for the southern relic. In the first 
case, the radio power is matched with the modest field strength of 
B le \i C ~ 0.7 /iG, while in the weaker southern shock the required 
magnetic field is « 219 fiG. Re-accelerated electrons can explain 
the observed radio power in both cases using fields in the range 
5reiic ~ 1 — 2/rG. The 7 -ray emission (lower panel) is the 

integrated hadronic emission in the downstream. Hence, the last 
bins on the left and on the right give the total emission from the two 
downstream regions and should be compared with the Fermi limit 
for A3667. When the contribution from both relics is summed up, 
the predicted hadronic emission is very close to the single-object 
limit from Fermi for this cluster. We find similar results for the 
relics in ZwCLJ2341 (see Tab. 2). In the next sections we will 
discuss the predicted magnetic fields and the 7 -ray emission from 
the full dataset and under different assumptions in our model. 


The full range of estimates from our fiducial model (Sec. 
m 3 is shown in Figure [4] where we plot: a) the distribution of 
predicted 7 -ray emission for the full sample (histograms with 
different colours for each model) and compare the mean emission 
of each run (thin lines) with the limits we derive from our stacking 
of Fermi exposures on these objects (Sec. 12.2.21 ). This stacking 
gives us the most robust testing of DSA, since it comes from the 
same set of objects simulated with our semi-analytical method. 
In the same figure, we also show for completeness the result of 
stacking only non-cool-core ( NCC) clusters in a larger sample of 
objects observed by Fermi ( iHuber et ali 20131) . which we converted 
into the [0.2 — 300] GeV energy range by assuming a 7 -ray 
spectral index of —2.6. The limit he re i s ~ 5 times lower than 
our stacking limit (see also lGriffin et al.lf2014l for a slightly lower 
limit), given the larger sample of objects (32) and the fact that 
this sample contains more nearby objects than our list of radio 
relics. This limit comes from a bigger population, yet comparing 
it to our simulated population is useful, under the hypothesis that 
two parent population of objects are dynamically similar. This 
is likely, because to the best of our knowledge all objects of our 
sample are NCC and all sho w evidence of very p erturbed dynam¬ 
ical st ates in X-ra y (e.g. |EdgeetaL 2003); Cavagnolo et alJ 
120091: iBagchi et al] 1201 ll: iBonafede et aD I 2 OI 2 L see also 
http://www.mergingclustercollaboration.org/merging-clusters.html). 
The magnetic fields at each relic required by our modelling of 
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the radio power, using EqQ] We also show for comparison the 
upper limits derived on the magnetic field at the location of relics 
in the Coma cluster feonafede et alJl2013l . based on the analysis 
of Faraday Rotat ion) and in the cluster CIZA 2242.2+5301 
l lvan Weeren et al] 120101 based on the analysis of the brightness 
profile across the relic), as well as the range of values inferred by 
iFinoeuenov et alj ( l2010l) for the relic north of A3667, based on the 
lack of Inverse Compton emission. 

In our model, both e = 0 and e = 5% runs predict a very 
high mean level of hadronic emission, well above the stacking of 
this dataset in the single injection case, and just below it in the case 
of e = 5% (but larger than the stacking of the full Fermi catalog). 
In both cases the mean emission is kept at a high value by 1/5 
of bright objects in the sample (A3667 and ZwCLJ2341), however 
the bulk of all remaining objects is also characterised by emission 
of the the order of the full stacking of non-cool-core clusters. On 
the other hand, the e = 1 % model predicts a mean emission be¬ 
low the the stacking of this dataset. However, the magnetic fields 
required in this case are very large. Here 6 out of 20 relics require 
7?reiic 10 pG where the radio emission is detected, which is at 
odds with the few available estimates of magnetic fields from ob¬ 
servations (gray arrows). In a general sense, explaining fields larger 
than a few ~ pG at the large radii where these relics are found 
is difficult based on both observational and theoretical facts (see 
Sec. l3.2.fl for a detailed discussion). 

To summarise, none of the scenarios we investigated with our semi- 
analytical method is able to simultaneously predict a level of 7 -ray 
emission compatible with Fermi and to use a reasonable magnetic 
field level in all objects. The reacceleration model assuming 1 % of 
CR energy to be reaccelerated by merger shocks survives the com¬ 
parison with Fermi, but makes use of very large magnetic fields in 
~ 1/3 of our objects. The models with single-injection or signifi¬ 
cant re-acceleration instead predict a mean emission for the sample 
which is in tension with stacking of Fermi observations for this 
sample, or for the larger sample of non-cool-core clusters, which 
very likely has the same characteristics of the double-relics one. 

The following sections will show how this problem is wors¬ 
ened as soon as all relevant assumptions made in our modeling are 
relaxed. 


3.2 Model uncertainties 

3.2.1 Magnetic field 

First, we investigate the role played by the magnetic field in relics. 
Explaining magnetic fields larger than a few /iG outside of cluster 
cores is very difficult for several reasons. 

In the only case in which a good volume coverage of the ICM 
is obtained through Farada y Rotation (i.e. in the Coma cluster, 
iBonafede et aljEoidl . 1 201 jl) . the inferred trend of magnetic field 
is B ~ n _ “ B , with q;b ~ 0.5 — 0.9, which implies that on average 
the field drops below 1/iG a t half of the viria l radius. This scaling is 
supported by simu lations JPolag et alJI 19991 : IBonafede et al.ll 20 ld : 
IVazza et alj2014ah and it implies that the magnetic field in the ICM 
is not dynamically important, i.e. /3 p i ~ 100 everywhere (where 
/3 p i = nk^T/P b and Pb is the magnetic pressure). Instead, a field 
of the order ^ lOpG at half of the virial radius or beyond implies 
/3 ^ 1, which is hard to justify theoretically. Indeed, the turbulence 
around the relic should be modest and dominated by compres¬ 
sive modes ; _and_cmi_only raise the magnetic field by a small fac¬ 
tor Jlanichino & Briiggenll201 AlSkillman et alj|2013l) . It has been 


suggested that CRs can cause magnetic field amplification via CR- 
driven turbulent amplification iBruggenll2013l) . but not to the ex¬ 
treme level required by our modelling. Moreover, all cosmologi¬ 
cal simulations predict some local amplification at shocks, but this 
is always smaller than the steep increase of the gas ther mal pres- 
sure, due to Rankine-Hugoniot jump conditions ( Dolag et al .1991 

D9r IVazza et al.ll2 


iBruggen et ai]|2005l: IXu et aljEool IVazza et al.ll2014ah . Based on 

simulations, the observed mass-radio Luminosity relation for dou¬ 
ble relics is better explained by assuming mag netic fields of the 
order of ~ 2 /iG in most objects dde Gasperin et alj|2014l ). 
Evidence of polarised radio emissions from a few radio relics (in¬ 
cluding a few contained in our sample here) exclude the pres¬ 
ence of large 10 pG fields distributed in scales below the ra¬ 
dio beam (a few ~ kpc), which would otherwise totally depo¬ 
larise the emission I van Weeren et al.ll 20 ld : IBonafede et alJl 20 lH . 
120121) . The highest (indirect) indication of magn etic fields in a gi- 
ant_ radio relic so far is of the order of ~ 7 pG ( Ivan Weeren et al J 
l20ldh . Moreover, all observations of Faraday Rotation outside of 
cluster c ores are consistent with magnetic fields of a few pG 


at most (Murgiaetal 


l20ldl: V; 


_ _ 20041 : iGuidetti et al J 120081 : IBonafede et alj 

et al J201C ). Finally, we notice that the acceleration ef¬ 


ficiencies usually assumed within PSA (e.g . iKang & Jonesll2007l : 
iHoeft & Briiggerj 120071 : iKang & R vul 120 1 3fi are based on the as¬ 
sumption of shocks running in a high /3 plasma. In the case of 
high magnetisation or relics, /3 -C 1, the physics of the intracluster 
plasma changes dramatically and the efficiencies from PSA are not 
applicable anymore. Hence, we test a scenario in which we cap the 
magnetic field at 33 re iic,max = 10 pG. In the many cases where the 
magnetic field inferred from Eq[T] would be larger than the 10 pG 
upper limit, we allow for the presence of reaccelerated electrons 
and iteratively increase e (by 0.1 % at each iteration) so that the ac¬ 
celeration efficiency is increased. We stop the iterations when the 
radio emission from Eq[T]matches the radio power (within a 10 % 
tolerance). We then assume this ratio to be constant in the down¬ 
stream region 0- Fig- 0 shows our results (see also Table 2). All 
magnetic fields are now more in line with the range of uncertainties 
given by the (scarce) observational data. In this case, the difference 
between all models is reduced, given that the energy ratio e had to 
be increased in several objects also when the starting model is the 
single-injection one. The model producing the largest emission is 
the one with e = 5%, but the difference with the two others is now 
limited to a few percent in the 7 -ray flux. In this case, all models are 
now at the level of the observed stacking for this dataset, and a fac¬ 
tor ~_2 abov e the full stacking of non-cool-core clusters in Fermi 
dlluber et alj|2013l ). We think that this set of runs gives the most 
stringent test to the PSA model because it includes the effect of 
CR re-acceleration to explain radio re lics for modest magnetic field 
dKang et alJl 2012 l : |Pinzke et al.ll2013l) . yet this idea fails when also 
the hadronic emission from accelerated CRs is taken into account. 
In the following we will discuss the remaining model uncertainties 
based on the ”Bcap” model. 


3.2.2 Upstream gas density and temperature 

Our assumption for the upstream gas density follows from the 
simplistic assumption of a /3-model profile along the direction 


5 In this case, the values of e = 0, e = 1% and e = 5% quoted in the 
labels actually refer to the initial assumed value for each cluster, while the 
final value depends on the iterations described here. However, in all cases 
the iterations stopped after reaching e ~ a few percent. 
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Table 2. Forecast of hadronic 7 -ray emission from the downstream region of our simulated clusters in the [0.2 — 300] GeV range. The first three columns 
show the prediction for the three assumed ratios e, and without imposing a maximum magnetic field at relics. The second three columns show our predictions 
by imposing a B ^ 10 fiG cap to the magnetic field, (see Sec. imi for details). For comparison, we show the upper limits from Fermi given in Table 1 for the 
same energy range. We mark with stars the objects for which the single-object comparison with Fermi data is problematic in several models. 


model 


object 


e = 0 e = 0.01 e = 0.05 e = 0, B c & P e = 0.01, B c &p e = 0.05, B ca , p observed 

lo gl0( e 7) lo Slo( e 7) lo glo( e 7) lo glC>( e 7) lo glo(e 7 ) l°glo( e 7) lo gl0 (^Fermi) 

[ph/(scm 2 )] [ph/(scm 2 )] [ph/(scm 2 )] [ph/(scm 2 )] [ph/(scm 2 )] [ph/(scm 2 )] [ph/(scm 2 )] 


A3376 

-10.07 

- 10.54 

-9.83 

-9.78 

-9.77 

-9.77 

-8.65 

A3 3 65 

-12.52 

- 12.12 

11.41 

-11.94 

-12.02 

-11.31 

-8.61 

A1240 

-10.67 

-11.60 

-10.90 

-10.84 

-10.83 

-10.84 

-9.21 

A2345 

-10.24 

-10.75 

-10.05 

-10.00 

-9.98 

-9.98 

-8.55 

RXCJ1314 

-10.69 

-10.94 

-10.24 

-10.18 

-10.16 

-10.16 

-9.06 

MACSJ1149 

-11.31 

-12.24 

-11.53 

-11.47 

-11.46 

-11.46 

-9.53 

MACSJ1752 

-10.14 

-10.45 

-10.86 

-10.80 

-10.75 

-10.79 

-8.92 

A3667* 

-8.53 

-9.28 

-9.28 

-9.23 

-9.21 

-9.20 

-8.59 

ZwCLJ2341* 

-8.80 

-9.81 

-9.11 

-9.05 

-9.04 

-9.04 

-9.41 

PLCKG287 

-10.75 

-11.44 

-10.74 

-10.68 

-10.67 

-10.67 

-8.81 


basic 

model 
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Figure 4. Left: distribution of predicted 7 -ray emission from our cluster sample, for different fiducial models (colored histograms). The thin vertical lines 
show the mean emission for the sample according to each model, and should be compared with the upper limits from the stacking of all non-cool-core clusters 
observed by Fermi, or by the stacking limited to our sample of clusters with double relics (gray arrows). Right: distribution of magnetic fields (colored 
histograms) required for each relic to match the observed radio power using Eq. |T| The additional gray arrows give the range of values in (erred for the few 
observations of magnetic fields in real relics (see text for details), while the hatching marks the values of magnetic fields that we regard as physically unrealistic 
(see Sec l3.2.lV 


of propagation of merger shocks. However, clusters hosting dou¬ 
ble relics are known to be perturbed. Observations provide evi- 
denc e that radio relics are a ligned with the merger axis of clus¬ 
ters l lvan Weeren et al.ll201 lbl ) and cosmological simulations show 
that the gas density along the major axis of merging clusters is 
significantly higher than the average profile, up to ~ 20 — 30 % 
close to the virial radius dVazza et alJl 201 lbl : iKhedekar et all201~3l : 
IVazza et ~aTl 120 1 3l ). X-ray observations a lso suggest departures of 
this level in the outer parts of cluster s feckert et~aill 2 QI 2 L 1 201 3l ; 
lUrban et al.ll2014l : iMorandi et ap2015l) . 

We tested the impact of a systematically 20 % higher upstream 
gas density in all our clusters, producing the results given in Fig. 
[ 6 ] (left). The enhanced density exacerbates the problems with the 
7 -ray emission because this scales as oc n 2 . The localised presence 
of denser clumps along the major axis of relics can only make this 


problem worse. We conclude that the hadronic emission predicted 
by our baseline model probably underestimates the level of 7 -ray 
emission that DSA should produce in these objects. Moreover, our 
assumption on the upstream gas density can be relaxed by con¬ 
sidering that very likely before the heating by the crossing merger 
shocks the medium in front of the relic had a temperature ^ T 500 • 
As a very conservative case, we considered a pre-shock tempera¬ 
ture lower by a factor ~ 2 compared to T500, corresponding to a 
M ~ 2 shock. The right panel of Fig. [ 6 ] shows the result of this 
test. The predicted 7 -ray emission is somewhat reduced compared 
to the ”Bcap” model, most notably in the single-injection case be¬ 
cause the propagation history of the single shock is the only crucial 
parameter. In this case the mean emission is a factor ~ 2 below 
the Fermi limits. However, in this case even more relics in all mod¬ 
els require ~ 10/iG fields in more objects because a decrease in 
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maximum B=10/xG 
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Figure 5. Similar to Fig[4] but here we impose a maximum magnetic field of 10 //G at the location of relics, and allow CR re-acceleration in all models (See 
Sec. l3.2.Tl for details). 


the upstream temperature in Eq[T|must be balanced by an increased 
magnetic field to match the radio power. In a more realistic case, we 
expect that the two above effects are combined since large-scale in¬ 
fall pattern along the major axis of merger clusters push cold dense 
un-virialized material further into the virial radius of the main halo, 
and therefore the problems of our DSA modeling of radio relics 
should become even worse. 

3.2.3 Radial dependence of the Mach number 

Our assumption of a constant Mach number within the volume 
of interest follows from the assumption that the kinetic energy 
flux across shocks is conserved during their propagation. This 
implies that the dissipation of kinetic energy by these shocks is 
negligible, which is appropriate for the weak shocks considered 
here. Also the upstream medium is assumed to be isothermal (i.e. 
n n (r)(Mc s ) 3 S(r)/2 = constant which gives M oc [n u (r) • 
S(r )] -1 / 3 ss constant because n u (r) cx r ~ 2 outside of the clus¬ 
ter core in the /3-model). We also tested the possibility of a shallow 
radial dependence of the Mach number, M(r) = Mo(r/_R re ii c ) 1 / 2 
(where Mo is the Mach number estimated at the location of the 
relic), as ir l Vazza & Briiggenl ( l20l4) . This was derived from the ob¬ 
served radial dependence of the radio spectral index of relics with 
radius ( Ivan Weeren et alj|2009t) . However, when only giant radio 
relics are considered, this tren d is not significant dBonafede et alj 
12011 ide Gasperin et alJl20l4h . To the best of our knowledge the 
average radial trend of Mach n umber in clusters wa s discussed 
only b y lvazza et al] d2009L [ 201 oh and more recently bv lHong et alj 
d20 1 4]) . All these works confirm a very shallow functional depen¬ 
dence with radius of the average Mach number of shocks, typically 
going from M ~ 1.5 — 2 in the centre to M ~ 3 in the clus¬ 
ter periphery. This trend is consistent with M oc Mo(r/rf) 1 ^ 2 
(with Mo ~ 2). In Fig. [7] we show the results if we impose a 
M oc r 1 ' 2 instead of a constant Mach number. The predicted 7 - 
ray emission downstream of relics is significantly reduced only in 
the single injection case (e = 0 ), but the average emission of the 
sample still remains larger than both stacking limits. We conclude 


that the radial trend of Mach number with radius, as long as this is 
shallow as suggested by simulations, is not a crucial point. In the 
shock re-acceleration case, we also tested a scenario in which the 
re-acceleration is done by a M = 3 shock instead of the M = 2 as 
in our baseline model (Fig. [7] right). In this case the requirement on 
the magnetic field are lowered and the number of relics for which 
we require B re u c ^ 10 pG is limited to one object (not shown). 
However, the predicted level of hadronic emission is now much in¬ 
creased and also the e = 1% run hits the Fermi stacking limits for 
this sample of objects, while the e = 5% now predicts an average 
emission which is ~ 4 — 5 times larger than this. 

3.2.4 Mach number from the radio spectrum 

The estimate of the Mach number at the position of relics is crucial 
in our modeling, as it determines the level of CR acceleration in the 
DSA scenario we are testing. However, several effects can make the 
Mach number we derive from the radio spectrum in Sec. [2T] un¬ 
certain. In several observations the Mach number derived from the 
radio spectrum is_foundto be higher t han th e one estimated through 
X-ray analysis fogrean et alJl2013l : lAkamatsu & Kawaharall2013h . 
The surface of complex shocks is described by a range of values 
of Mach numbers rather than by a single value ( iSkillman et alj 
120131 ) and the radio emission will be dominated by e lectrons p rob- 
ing larger Mach numbers compared to the mean l lHong et alj2014l) . 
Additionally, radio observations with only a few beams across the 
relic can only produce integrated radio spectra, which are expected 
to be by « 0.5 steeper than the injection spectrum. A blend of 
several populations of electrons seen in projection can yield spec¬ 
tra with time-dependent biases, as recently discussed by [Kang 
( i2014l) . To assess this effect, we run a set of Monte Carlo meth¬ 
ods and extract uniform random deviates within M ± AM, where 
AM ^ 0.5 M. For each relic we randomly extracted 200 values 
of AM and compute the downstream 7 -ray emission in all cases. 
Figure[ 8 ]shows the results for the single injection and reacceleration 
models (colored histograms). In this case the simulated mean emis¬ 
sion in the plot is the average of each set of 200 realisations (i.e. we 
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20 % larger upstream density 



50% lower upstream temperature 



Figure 6. Left panel: distribution of predicted 7 -ray emission from our cluster sample, similar to Fig. [4]but here assuming a 20 % higher upstream gas density 
(left) compared to the basic model. Right panel: as in the left panel but assuming a 50 % lower upstream gas temperature compared to the basic model (see 
text for details). In both cases we assume a capping of the magnetic field at 10 fiG. 


M scales with radius 



reacceleration from Mach = 3 



icr 13 icr 12 icr 11 icr 10 icr 9 icr 8 

7 —emission[ph/(s cm 2 )] 


Figure 7. Left: distribution of predicted 7 -ray emission from our cluster sample, similar to Fig. [5] but here considering a radial scaling of the Mach number 
downstream of relics, M(r) = Mq (r / /Pelic :) 1 f 2 ■ Right: same as in Fig. [4] but here assuming that the re-acceleration in the e = 1% and e = 5% cases is 
done by a M = 3 shock. 


first compute the mean emission for the cluster sample, for one ran¬ 
dom combination of extractions for AM, and then compute the av¬ 
erage emission and dispersion within the full dataset of 200 random 
realisation). The number of 200 realisation was chosen based on the 
fact that the errors in the mean emission do not change significantly 
for larger numbers of realisations. Compared to our fiducial model, 
a random variation in the assumed Mach number overall increases 
the level of hadronic 7 -ray emission, suggesting that on average 
our baseline model underestimates the total emission. The under¬ 
estimate is obviously more significant in the single-shock model 


(e = 0), where the mean emission is ~ 5 smaller than what we ob¬ 
tain as an average from the 200 random extractions. In the reccel- 
eration cases, the effect is smaller, and the fiducial model probably 
underestimates the hadronic emission by a factor ~ 1.5 — 2. For 
the same set of random extractions, the problem with requiring too 
large magnetic fields for a significant fraction of objects (~ 1/3 
in the e = 0 and e = 1 %) remains (not shown). We conclude 
that realistic uncertainties in the Mach number do not change the 
robustness of our results. 
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random Mach number 



7-emission[ph/(s cm 2 )] 



10-13 10 -i2 10 -11 10~ 10 10~ 9 10~ B 1CT 7 
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Figure 8. Distribution of predicted 7 -ray emission from our cluster sample 
including a random deviate from the Mach number derived from the radio, 
in the range M = M ra( j [ 0 =b 0.5M raf ji o (see Sec. Hl4l for details). We ex¬ 
tracted 200 random deviates for each object and computed the downstream 
hadronic emission for the three reacceleration models. Differently from the 
previous Figures, in this case the vertical lines show the mean emission 
from the 200 stacked samples (i.e. we computed the mean emission within 
the sample, one time for each random extraction, and then computed the 
mean emission over the 200 realisations). 


3.2.5 Viewing angle 

So far we assumed that all relics trace shocks which propagate ex¬ 
actly in a plane perpendicular to the line of sight. Numerical simu¬ 
lations of relics support this scenario and limit the inclination along 
the lin e of sight of the propagation plane down t o ^ 10 — 20 de¬ 
grees l lvan Weeren et alfcoi lallKanp et all2012l) . In general, sim¬ 
ulated radio relics assuming DSA resemble the observational prop- 
erties of most relics only when they lie cl ose to the plane of the sky 
dVazza et al.ll2012bl : ISkillman et all201~3h . However, very small in¬ 
clinations can be present and we checked if the inclusion of small 
(| Aw| 7 30 degrees) along the line of sight can alter the picture in 
any significant way. Similar to the previous test, we randomly ex¬ 
tracted 200 values of Aw for a uniform distribution in the ”Bcap” 
model, and accordingly recalculated the total volume spanned by 
shocks (which can only become bigger compared to the Aw = 0 
case), and computed the average value of the 200 realisations of 
cluster stackings. Fig. [9] shows the results of this test. The aver¬ 
age 7 -ray emission from all realisations is of the order the fiducial 
model (Aw = 0) and at the level of the Fermi s tacking for thi s 
sample, and larger than the stacking by lAckermann et alj ]20l4) . 
The outcome in the distribution of magnetic fields at relics is even 
worse than in the fiducial case, because in the case of large an¬ 
gles along the line of sight the relics are located further out, where 
the gas density is lower than in the Aw = 0 case, and the mag¬ 
netic field must increase dramatically to match the radio power. As 
for all previous tests, we conclude that the presence of small but 
unavoidable projection effects has a small effect. However, on av¬ 
erage these projection effects should yield an even larger hadronic 
emission from our dataset if DSA is at work. 


Figure 9. Similar to Fig[ 8 ] but here by extracting 200 random values for the 
angle in the plane of the sky for each relic, |<5w| 7 30 degrees. 


3.2.6 Uncertainties in Cosmic Ray physics 


The efficiencies that we have tested produce the lowest amount 
of CR-acceleration dKang & Rvul 120131 ). However, previously 
suggested funct ions for the t)(M) acceleration efficiency (e.g. 
iKang & Jonesll2007l) give a larger injection of CRs from M 7[ 5 
shocks, and can only make the problem with Fermi limits worse. 
Other uncertainties in the physics of CRs after their injection are 
briefly discussed here. Outside of cluster cores, the CRs are only 
weakly subject to hadronic and Coulomb losses, owing to the low 
gas density. For the sake of our analysis, it does not make any differ¬ 
ence if they diffuse in the cluster volume at constant radius (since 
their contribution to the 7 -ray emission only depends on the gas 
density and not on the exact location in the cluster atmosphere). 
Only diffusion in the vertical sense can change our estimate. How¬ 
ever, this cannot be a big effect since CRs are thought to be frozen 
into tangled magnetic fields, and in this case their spatial diffusion 
is slow (i.e. r ~ 2 ■ 10 8 yr (R/Mpc ) 2 (E/GeV ) -1 ^ 3 for a con¬ 
stant B = 1 pG m agneti c field and assuming Bohrn diffusion, 
e.g. [Berezinsky et all 19971) . More recently, it has been suggested 
that the fast (^streaming ^ va, where va is the Alfven velocity) 
streaming of CRs can progressively deplete the downstream re¬ 
gion of the shock and reduce radio and hadronic emission (e.g. 
EnBlin et alj |201ll) . offering a way to reconcile hadronic models 
for radio halos with the observed bimodality in the distribution of 
diffuse radio emission in clusters dBrunetti et al]|2007t>. However , 
the validity of this scenario is controversial (e.g. Donnertll2013l : 
Wiener et al 1 120131) and this mechanism has been suggested to be 
maximally efficient in relaxed clusters. Instead, our sample of clus¬ 
ters with double relics is all made by objects with a very unrelaxed 
X-ray morphology. In the case of_turbulent intracluster medium, 
the detailed calculation bv lwiener et al.l ( f2013t) shows that even fast 
streaming can also rapidly diminishes the 7 -ray luminosities in the 
E = 300 — 1000 GeV energies probed by imaging air Cerenkov 
telescopes (MAGIC, HESS, VERITAS ), but not in the lower energies 
probed by Fermi. Therefore, the latter is a more robust probe of the 
CR injection history, and our modeling here is well justified. Fi- 
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nally, our work neglects the (re)acceleration of CRs by other mech¬ 
anis ms, such a s turb ulent reacceleration in the downstream of relics 
(e.g. lBrunetti et all2004h . which should re-energise radio emitting 
electrons as well CR-protons. Including this effect in our simple 
modeling is beyond the goal of this work, however its net effect 
can only be that of further increasing the mean emission we predict 
here. 


3.3 Other observational proxies: inverse Compton emission 


Relativistic electrons accelerated by shocks can also emit in the 
hard-X ray band through Inv erse Compton (IC) emission (Rephaeli 
Il979l : lsarazin & Lieull 19981) . In principle, this can offer a comple¬ 
mentary way of testing our models without having to make assump¬ 
tions for the magnetic field. Hence, we have computed the IC emis¬ 
sion from each object in our sample, under th e ass umption of sta¬ 
tionary shock acceleration following[Sarazi rjj (t 1999lf: 

eic«0.17^^( 7 ^5-10 3 ), (8) 

where At is given by the shock crossing time for each radial shell, 
EcR,e is the energy of CR-electrons injected by the shock, and our 
integration is limited to the cooling region close to each shock. This 
is computed using Eq. (12) in lKang et al.l ( l2012h : 


(cool ~ 890 kpc 


Vs 


B 1 / 2 


(1 + z) —1/2 , (9) 


10 3 km/s B 2 + B 2 MB 1 GHz 


where z is the redshift. Here we show our prediction for the last 
investigated case, where we cap the magnetic field at B = 10 fiG 
and allow for the inclusion of CRs re-acceleration when necessary 
to match the observed radio power (Sec. refsubsubsec:Bfield). 

Figure[l0]gives our predicted emission for the single injection 
case (blue) and for the extreme re-acceleration case (e = 0.05), 
to emphasize the weak dependence of the predicted IC emission 
on the assumed re-acceleration model. Table 3 gives the predicted 
flux in IC emission from the downstream of all double relics in 
the [20 — 100]keV range, and the assumed radiative lengths in 
the downstream region. In all cases the predicted emission lies 
below the detection threshold by the hard-X ray satellite NUS- 
TA/ 0 , which has been estimated to be of the order of a few 
~ 10 -12 erg/(s cm 2 ) in the case of th e recent observations of 
the C oma cluster dGastaldello et al.ll2014l) and of the Bullet cluster 
dWik et all2014l ). However, a significantly lower sensitivity might 
be reached in the case of peripheral relics, given that in the latter 
clusters the contamination from the hot thermal gas in the hard-X 
ray range hampers the detection of the IC signal. This might be 
the case for the most powerful targets in our sample, represent by 
A3376 (both relics) and by the most powerful relic in ZwCLJ2341. 
In these cases, the large distance from the centre of host clusters 
(~ 1.2 — 1.3 Mpc) might indeed offer a better chance of detection 
of the IC signal. In the next years, the Astro-H satellite should be 
able to probe the i nverse Compton em ission in the same clusters 
dAwaki et alJl2014l : lBartels et all 201 5l) . 


4 DISCUSSION AND CONCLUSIONS 

Most of the evidence from X-ray and radio observations suggests 
a link between radio relics and merger sho cks: merger axes 
of clusters and relic orientations correlate dvan Weeren et all 

6 http:www.nustar.caltech.edu 



Figure 10. Forecast of Inverse Compton emission from our simulated relics, 
in the extreme cases of e = 0 (blue) and e = 0.05 (red). 


1201 lb ), the power of r adio r elics scales with X-ray lumi¬ 
nosities _|BonafedeetalJ |2012l : iFeretti et al.l 12012l) and mass 
(de Gasperi n et al. 2014 ) of the host cluster. Cosmological sim¬ 
ulations produce emission patterns consistent with observed 


across shock waves (e.g. 

Hoeft et al. 2008i; Pfrommer 

2008 

Battaglia et alj 2009|; |skillman et al. 

201 ll: VazzaetalJ 

2012 b 

Nuza et al. 2012: Skillman et al. 2013 

). 


Still, a number of recent observations have revealed 
som e 7 >pen issues, including uncertain merger scenarios (e.g. 
QgreanetalJ 201 31). departures from power-law spectra (e.g. 


Stroe et al 20141 : iTrasatti et afl 120141) . missi ng associations be¬ 


tween radio emissi on and X-ray maps (e.g. iRussell et al.l 1201 ll: 
Qgrean et afl 120141), effici encies problems (e.g. iMacario et all 


2011 ; Bonafede et al J|20 1 211 . inconsistencies between Mach num¬ 


bers derived from X-ray and radio observations (e.g. Kale et all 


120121 : lAkamatsu & Kawaharall20131 : lOgrean & Bragged 201 3l) and 

apparent connections to radio galaxiesi lBonafede et al.ll2014 ). 


In this work, we used a simple semi-analytical model of ex¬ 
panding merger shocks in clusters to reconstruct the propagation 
history of shocks leading. We used a spherically symmetric model 
and assumed that cosmic ray protons are trapped in the intracluster 
medium on all relevant timescales. A range of realistic scenarios for 
the acceleration of relativistic electrons and protons via DSA, vary¬ 
ing the upstream gas conditions, the shock parameters and the bud¬ 
get of pre-existing cosmic rays, gives very similar results. In all re¬ 
alistic scenarios, a significant fraction of our objects (~ 1/2 — 1/3) 
has difficulties in matching at the same time the observed radio 
emission and the constraints imposed by the Fermi limits, unless 
the magnetic field in all problematic objects is much larger than 
what usually considered realistic (3> 10 [i G). The scenario in 
which radio emitting electrons comes from the re-acceler ation of 
pre-existing electrons ( Ka ng et al.l2012l : |Pinzke et al.l20l 3ll can al¬ 
leviate the tension with Fermi if the pre-existing electrons are not 
the result of previous injection by shocks as we investigated here, 
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Table 3. Forecast of IC emission from the downstream region of relics in our simulated clusters, for [20 — 100] keV. The predictions are here only given for 
the e = 0 using our model with a capping of the magnetic field at 10 fiG, see Sec. non Tor details), as all others yield extremely similar results. 


object Ml M 2 log 10 (eic,i) l°g 10 ( e ic,2) Zcool.i (cool, 2 B 2 

[erg/(s cm 2 )] [erg/(scm 2 )] [kpc] [kpc] [/xG] [/xG] 


-13.98 -13.19 


A3376 

3.3 

3.3 

A3365 

2.1 

1.8 

A1240 

3.3 

2.8 

A2345 

2.8 

2.2 

RXCJ1314 

2.4 

2.4 

MACS J1149 

3.3 

2.4 

MACSJ1752 

3.2 

4.0 

A3667 

1.7 

3.8 

ZwCLJ2341 

1.8 

3.2 

PLCKG287 

2.9 

2.2 


-16.78 

-18.78 

-15.29 

-15.24 

-14.41 

-15.40 

-15.11 

-15.28 

-15.54 

-16.29 

-14.88 

-14.92 

-14.53 

-14.12 

-16.96 

-13.16 

-14.64 

-16.05 


273.6 

191.3 

2.25 

0.54 

87.6 

48.6 

5.41 

8.69 

100.7 

96.5 

0.99 

0.89 

107.2 

131.1 

0.60 

0.97 

86.4 

111.2 

0.38 

0.65 

101.7 

96.4 

1.75 

1.46 

122.2 

89.1 

3.21 

5.92 

164.4 

159.7 

1.56 

1.32 

185.7 

97.7 

2.15 

0.36 

125.4 

128.4 

1.40 

3.89 


but are instead released by mechanisms that mostly inject leptons 
(e. g. leptonic-do min ated jets from AGN), as already discussed in 
IVazza & Briiggenl (120141) . 

Based on our semi-analytical model, the standard DSA 
scenario with thermal leakage that predicts that Ecn e -C -Ecitp 
cannot simultaneously explain radio relics and produce less 
7 -radiation than the upper limits from Fermi , unless unrealistically 
large magnetic fields are assumed at the position of relics (e.g. 
Brelic ^ 10 — 100 /xG). This result is very robust, at least in the 
statistical sense, against all investigated variations of our fiducial 
parameters for the modeling of the shock acceleration of CRs. 
Additional effects that go beyond our idealized modeling of cluster 
mergers, e.g. a clumpy ICM, a succession of mergers and the 
additional acceleration of CRs by AGN, supernovae, turbulence or 
reconnection exacerbate this discrepancy. 

Despite its obvious degree of simplification, a semi-analytical 
method is useful to tackle the case of double relic systems. In 
these systems it is reasonable to assume that most of the energet¬ 
ics is related to the observed pair of giant merger shocks. Their 
shape and location is rather regular and symmetric with respect to 
the cluster centre, suggesting that one can make reasonable esti¬ 
mates for their propagation history. This setup allows us to run very 
fast testing of different possible acceleration scenarios, and as we 
showed in our various tests it generally gives a lower limit on the 
expected 7 -ray emission. This method is meant to be complemen¬ 
tary to fully cosmological numerical simulations where the effects 
of multiple shocks, particle advection and cooling, as well as inputs 
from galaxy formation and other mechanisms can be taken into ac - 
count at run-time (e.g. IPfrommer et alj|2007l : IVazza et all 1 201 2 al) . 
However, a thorough exploration of models is computationally de¬ 
manding because of the required high resolution and the complex¬ 
ity of the numerics. Also the agreement between different numeri¬ 
cal techniques on this topic is still unsatisfactory (see discussion in 
I Vazza et aOl201 lal) . 

Another way of illustrating our result is found by rescaling 
the efficiency for proton acceleration, such that the upper 

limits from Fermi are not violated. In the case without pre-existing 
CRs, this is a simple exercise as we only need to rescale r/(M) 
for each relic separately, and compute the average of the efficien¬ 
cies for each bin of the Mach number (here we chose a bin size 
of AM = 0.6 to achieve a reasonable sampling of the sparse 
distribution of Mach numbers in the dataset). Here, we keep the 
magnetic field fixed_at B =_2uG as suggested by recent obser¬ 
vations dde Gasperin et al .1120141) . The result is shown in Fig. EU 
where we show the maximally allowed acceleration efficiency for 



Figure 11. Acceleration efficiency of CR-protons (green, rescaled by a fac¬ 
tor X10 down) and CR-electrons (blue), and electron to proton acceleration 
ratio (red) allowed by our combined radio and 7-ray comparison with ob¬ 
servations. In this case, we assumed a fixed magnetic field of B = 2jiG for 
all relics. 


CR-electrons, protons, as well as K e / p as a function of M. This 
relation results from a somewhat coarse simplification of the prob¬ 
lem but it is a rough estimate of the acceleration efficiencies in 
weak ICM shocks. 

For shocks with M ^ 2, the flux ratio of injected elec¬ 
trons is larger than that in protons, K e / p ~ 1 — 100 , at odds 
with standard DSA (even including re-accelerated electrons). For 
M > 2.5 the acceleration efficiency of protons can become signif¬ 
icant (~ 10~ 3 — 10 -2 ) while the acceleration efficiency of elec¬ 
trons flattens and K e / p ~ 10~ 2 . The functional shape (rf the accel¬ 
erat ion efficiency for protons is consistent with the iKang & Rvd 
120131) model, but the absolute normalisation is lower by a factor 
~ 10 - 100 . 

A possible solution has been suggested bv lKang et alJ J2014I) . 
who assumed that electrons and protons follow a ^-distribution 
near the shock transition. A xc-distribution is characterized by a 
power-law rather than by an exponential cutoff at high energies, 
thus ensuring a more efficient injection of high-energy particles 
into the DSA cycle. This distribution is motivated by spacecraft 
measurements of the solar win d as well as by observations of 
HII and planetary nebulae (e.g. Lazar et all 12012i and references 
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therein). iKang et all 120141 ) explored the application of the k- 
distribution to M ^ 2 shocks in the ICM, and concluded that the 
distribution can have a different high-energy tail as a function of 
the shock obliquity and of the plasma parameters. In the ICM, 
the distribution might be more extended towards high energies 
for electrons than for protons thus justifying a higher acceleration 
efficiency for electrons than for protons. However, in order to 
explain the origin of these wider distributions, one must resort to 
detailed micro-physical simulations of collisionless shocks. 

The most promising explanation for the non-observation of 
7 -rays has been suggested bv lGuo et all d2014l) who studied the ac¬ 
celeration of electrons with particle-in-cell (PIC) simulations under 
conditions relevant to merger shocks. They showed that M ^ 3 
shocks can be efficient accelerators of electrons in a Fermi-like 
process, where electrons gain energy via shock drift acceleration 
(SDA). The electron gain energy from the motion of electric field 
and scatter off oblique magnetic waves that are self-generated via 
the firehose instability. They found that this mechanism can work 
for high plasma betas and for nearly all magnetic field obliquities. 
However, these simulations have been performed in 2D, and could 
not follow the acceleration of electrons beyond a supra-thermal 
energy because of computing limitations. At the same time , hy- 
brid simulations of proton acceleration by ICaprioli & Spitkovskvl 
d2014l) have shown that the acceleration efficiency is a strong func¬ 
tion of the obliquity angle. If indeed the magnetic field in radio 
relics is predominantly perpendicular to the shock normal, as found 
e.g. in the relic in the cluster CIZA 2242.2+5301, then the predic¬ 
tion is that the acceleration efficiency of protons is strongly sup¬ 
pressed, thus explaining the non-detection of hadronic emission. It 
remains to be seen if the results of these simulations hold in 3D, 
with realistic mass ratios between electrons and protons and cou¬ 
pled to a large scale MHD flow. It is also not clear whether the 
magnetic field is quasi-perpendicular in all the relics of this sample 
and how the alignment of the magnetic fields with the shock sur¬ 
face observed on large scales can be scaled down to scales of the 
ion gyro radius. 
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